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Abstract. The production of n , rr - , K , K~, p, and p at mid- rapidity has been measured in proton-proton 
collisions at ^/s = 900 GeV with the ALICE detector. Particle identification is performed using the specific 
energy loss in the inner tracking silicon detector and the time projection chamber. In addition, time-of-fiight 
information is used to identify hadrons at higher momenta. Finally, the distinctive kink topology of the 
weak decay of charged kaons is used for an alternative measurement of the kaon transverse momentum (p%) 
spectra. Since these various particle identification tools give the best separation capabilities over different 
momentum ranges, the results are combined to extract spectra from p t = 100 MeV/c to 2.5 GeV/c. The 
measured spectra are further compared with QCD-inspired models which yield a poor description. The 
total yields and the mean p t are compared with previous measurements, and the trends as a function of 
collision energy are discussed. 



1 Introduction 



In pp collisions at ultra-relativistic energies the bulk of the 
particles produced at mid-rapidity have transverse mo- 
menta, pt, below 1 GeV/c. Their production is not calcu- 
lable from first principles via perturbative Quantum Chro- 
modynamics, and is not well modelled at lower collision 
energies. This low p t particle production, and species com- 
position, must therefore be measured, providing crucial 
input for the modelling of hadronic interactions and the 
hadronization process. It is important to study the bulk 
production of particles as a function of both p t and parti- 
cle species. With the advent of pp collisions at the Large 
Hadron Collider (LHC) at CERN a new energy regime is 
being explored, where particle production from hard in- 
teractions which are predominantly gluonic in nature, is 
expected to play an increasing role. Such data will pro- 
vide extra constraints on the modelling of fragmentation 
functions. The data will also serve as a reference for the 
heavy-ion measurements. 

The ALICE detector [TJ[2] is designed to perform mea- 
surements in the high-multiplicity environment expected 
in central lead-lead collisions at ^/snn = 5.5 TeV at the 
LHC and to identify particles over a wide range of mo- 



menta. As such, it is ideally suited to perform these mea- 
surements also in pp collisions. 

This paper presents the transverse momentum spectra 
and yields of identified particles at mid-rapidity from the 
first pp collisions collected in the autumn of 2009, during 
the commissioning of the LHC, at -Js = 900 GeV. The 
evolution of particle production in pp collisions with colli- 
sion energy is studied by comparing to data from previous 
experiments. 

We report 7r + , 7r~, K + , K - , p, and p distributions, 
identified via several independent techniques utilizing spe- 
cific energy loss, AE /Ax, information from the Inner Track- 
ing System (ITS) and the Time Projection Chamber (TPC), 
and velocity measurements in the Time-Of-Flight array 
(TOF) . The combination of these methods provides par- 
ticle identification over the transverse momentum range 
0.1 GeV/c < pt< 2.5 GeV/c. Charged kaons, identified via 
kink topology of their weak decays in the TPC, provide a 
complementary measurement over a similar p t range. All 
reported particle yields are for primary particles, namely 
those directly produced in the collision including the prod- 
ucts of strong and electromagnetic decays but excluding 
weak decays of strange particles. 

The paper is organized as follows: In Section 2, the AL- 
ICE detectors relevant for these studies, the experimental 
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49 conditions, and the corresponding analysis techniques are 

50 described. Details of the event and particle selection are 

51 presented. In Section 3, the 7r + , 7r~, K + , K~, p, and p in- 

52 elusive spectra and yields, obtained by combining the var- 

53 ious techniques described in Section 2, are presented. The 

54 results are compared with calculations from QCD-inspired 

55 models and the pt-dependence of ratios of particle yields, 

56 e.g. K/7r and p/71", are discussed. Comparisons with data 

57 from other experiments at different y^s are made and the 

58 evolution of the ratio of strange to non-strange hadrons 

59 with collision energy is discussed. Finally, in Section 4 the 

60 results are summarized. 



61 2 Experimental setup and data analysis 

62 2.1 The ALICE detector 

63 The ALICE detector and its expected performance are 

64 described in detail in [TjfJ]. For the analyses described 

65 in this paper the following detectors are used: the ITS, 

66 the TPC and the TOF detector. These detectors are posi- 

67 tioned in a solenoidal magnetic held of B = 0.5 T and have 

68 a common pseudo-rapidity coverage of —0.9 < r\ < 0.9. 

69 Two forward scintillator hodoscopes (VZERO) are used 

70 for triggering purposes. They are placed on either side of 

71 the interaction region, covering regions 2.8 < 77 < 5.1 and 

72 -3.7 < r? < -1.7. 

73 2.1.1 The Inner Tracking System 

74 The ITS is the closest of the central barrel detectors to the 

75 beam axis. It is composed of six cylindrical layers of sili- 

76 con detectors. The two innermost layers are equipped with 

77 pixel detectors (SPD), followed by two layers of drift de- 

78 tectors (SDD) and two layers of double-sided silicon strip 

79 detectors (SSD). The innermost layer is at 3.9 cm from 
so the beam axis, while the outer layer is at 43.0 cm. 

si The ITS provides high-resolution space points that 

82 allow the extension of tracks reconstructed in the TPC 

83 towards the interaction vertex, thus improving momen- 

84 turn and angular resolution. The four layers equipped with 

85 SDD and SSD also provide a measurement of the specific 
se energy loss dE/dx. The SPD yields an on-line measure 
87 of the multiplicity by counting the number of chips that 
as have one or more hits (fast-OR), which is included in the 

89 minimum-bias trigger logic [3]|3] . The ITS is also used as a 

90 stand-alone tracker to reconstruct charged particles with 

91 momenta below 200 MeV/c that are deflected or decay 

92 before reaching the TPC, and to recover tracks crossing 

93 dead regions of the TPC. A detailed description of the 

94 three sub-systems can be found in [3J. The dE/dx mea- 

95 surement in the SDD and SSD has been calibrated using 

96 cosmic ray data and pp events [S]. The 2198 ITS modules 

97 have been aligned using survey information, cosmic-ray 

98 tracks and pp data with the methods described in [5] . The 

99 fraction of active modules per layer in the present setup 

100 is around 80% in the SPD and 90% - 95% both in SDD 

101 and SSD. 



2.1.2 The Time Projection Chamber 102 

The TPC is the main tracking device. It is a large volume, 103 

high granularity, cylindrical detector with an outer radius 104 

of 2.78 m and a length of 5.1 m. The active volume extends 105 

from 0.85 m to 2.47 m in radius. It covers 2n in azimuth 106 

and |?7| < 0.9 in polar angle for the full radial track length. 107 

Accepting one third of the full radial track length extends ios 

the range to \i]\ < 1.5. The 90 m 3 drift volume is filled with 109 

a Ne (85.7%), C0 2 (9.5%), and N 2 (4.8%) gas mixture, no 

A high voltage central membrane splits the drift region in m 

two halves, resulting in a maximal drift time of 94 /is. Each 112 

of the two read-out planes is composed of 18 inner and 18 113 

outer chambers with a total of 159 pad rows, resulting m 

in a total of 557 568 pads which are read out separately. 115 

The position resolution in r<j> direction varies from 1100 u« 

/jm to 800 /jm when going from the inner to the outer 117 

radius. Along the beam axis (z, also the drift direction) us 

the resolution ranges between 1250 /im and 1100 /xm. A 119 

maximum of 159 clusters can be measured along a track 120 

in the TPC. For a detailed description see [7]. 121 

2.1.3 The Time-Of-Flight Detector 122 

The TOF detector consists of 18 azimuthal sectors, each 123 

containing 91 Multi-gap Resistive Plate Chambers (MR- 124 

PCs) distributed in five gas-tight modules. It is positioned 125 

at 370-399 cm from the beam axis. The region 260° < <fi m 

< 320° at 77 ~ is not covered in order to minimize the 127 

material in front of the Photon Spectrometer, which is not 12s 

used in this analysis. The MRPC detectors are installed 129 

with a projective geometry along the beam direction, min- 130 

imizing the variation of the flight path of particles across 131 

the sensitive area of the detector. Each MRPC is seg- 132 

mented into 96 read-out pads (2.5 x 3.5 cm 2 size), result- 133 

ing in a total of 152928 channels. Test beam results demon- 134 

strated that the intrinsic time resolution of the detector is 135 

better than 50 ps, dominated by electronic effects and the 136 

time resolution of the time-to-digital converters [8]. Re- 137 

suits from the TOF commissioning with cosmic rays are 13s 

described in references [9"HTl"]. In the present setup. 9.6% 139 

of the readout channels were inactive due to failures in 140 

the high- or low-voltage systems or in the readout elec- 141 

tronics. The fraction of noisy channels, identified during m 

data taking by online monitoring and excluded from the 143 

subsequent reconstruction, was below 0.1%. 144 

2.2 Event selection and normalization 145 

The data presented in this paper were collected during «6 

the commissioning of the LHC at CERN in the autumn of «? 

2009, with pp collisions at y/s = 900 GeV. The collider ua 

was run with four bunches per beam, resulting in two 149 

bunch crossings per beam circulation period (89 /is) at 150 

the ALICE interaction point. The remaining two bunches 151 

per beam were not collided at ALICE, and served to es- 152 

timate the contribution of beam-gas interactions. The av- 153 

erage event rate was a few Hz, so the fraction of pile- up 154 

events was negligible. 155 
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we The analysis is based on a sample of ~ 300k inelastic 

157 pp collisions. The online trigger selection requires a signal 

158 in either of the VZERO counters or at least one hit in ci- 

159 ther of the SPD layers. The selection was improved offline 
wo with recomputed trigger input quantities using the time 
lei average over all VZERO hits and a suppression of noisy 

162 channels. The contamination from beam-induced back- 

163 ground is rejected offline using the timing information of 

164 the VZERO and by cutting on the correlation between the 

165 number of clusters and track segments (tracklets) in the 

166 SPD detector [T2"1IT3"] . Selected events are further required 

167 to contain a reconstructed primary vertex. The vertex re- 

168 construction efficiency calculated via Monte-Carlo simula- 

169 tions is 96.5% for events with one reconstructed track and 

170 approaches unity for events with more than two tracks. 

171 The results presented in this paper are normalized to 

172 inelastic pp collisions, employing the strategy described in 

173 [T2, 13 . In order to reduce the extrapolation and thus the 

174 systematic uncertainty on the normalization, the sample 

175 of selected events used for normalization includes triggered 

176 events without reconstructed tracks or vertices. Those ev- 

177 ents still contain a small contamination from very low mul- 

178 tiplicity beam-induced background or accidentals from the 

179 trigger, which are not rejected by the selections described 

180 above. This contamination is of the order of 4% and is 
lsi subtracted using the control triggers. From the analysis 

182 of empty bunch events the random contribution from cos- 

183 mic rays is found to be negligible. The number of selected 

184 events is then converted to the number of inelastic col- 

185 lisions after correcting for the trigger efficiency, which is 

186 determined from the Monte-Carlo simulation, scaling the 

187 cross section for diffractive processes to the measurements 

188 of UA5 [H] . The subtraction of beam-gas events and the 

189 efficiency correction partially compensate each other: the 

190 overall correction factor is about 5% with a systematic un- 

191 certainty of about 2%, coming mainly from the uncertain- 

192 ties in the modelling of diffraction in the event generators. 

193 In order to compare to previous experimental results, 

194 which are only published for the non-single-diffractive 

195 (NSD) class, in Section^ we scale our spectra for the mea- 

196 sured ratio dN ch /dr]\ NSD / dN ch /d-q\ INEL ~ 1.185 [H]. 

197 PYTHIA and PHOJET simulations indicate that the p t - 

198 dependence of the ratio of spectra for NSD and inelastic 

199 collisions is less than 5% in the reported range. Particle 

200 ratios are found to be insensitive to the conversion from 

201 inelastic to non-single-diffractive events. 

202 2.3 Track selection 

203 The identified particle spectra were measured indepen- 

204 dently with the ITS, TPC and TOF, and combined in 

205 the final stage of the analysis. The rapidity range \y\ < 

206 0.5 was used for all analyses except for the kink analysis 

207 (\y\ < 0.7). 

208 For the TPC and TOF analyses, tracks reconstructed 

209 in the TPC are used. The TPC has full acceptance for 

210 tracks with \r}\ < 0.9. However, shorter tracks at higher 

211 rj can still be used for physics analysis, in particular pro- 

212 tons with a transverse momentum of p t = 400 MeV/c and 



| y | = 0.5 which correspond to |ry| = 1.1. To ensure high 213 

tracking efficiency and dE/dx- resolution, while keeping 214 

the contamination from secondaries and fakes low, tracks 215 

are required to have at least 80 clusters, and a \ 2 01 the 216 

momentum fit that is smaller than 4 per cluster. Since each 217 

cluster in the TPC provides two degrees of freedom and 21s 

the number of parameters of the track fit is much smaller 219 

than the number of clusters, the \ 2 cv ^ is approximately 220 

2 per degree of freedom. In addition, at least two clusters 221 

in the ITS must be associated to the track, out of which 222 

at least one is from the SPD. Tracks are further rejected 223 

based on their distance-of-closest approach (DCA) to the 224 

reconstructed event vertex. The cut is implemented as a 225 

function of pt to correspond to about seven (five) stan- 226 

dard deviations in the transverse (longitudinal) coordi- 227 

nate, taking into account the pt-dependence of the impact 22s 

parameter resolution. These selection criteria are tuned to 229 

select primary charged particles with high efficiency while 230 

minimizing the contributions from weak decays, conver- 231 

sions and secondary hadronic interactions in the detector 232 

material. The DCA resolution in the data is found to be 233 

in good agreement with the Monte-Carlo simulations that 234 

are used for efficiency corrections (see next Section). 235 

Tracks reconstructed in the TPC are extrapolated to 236 

the sensitive layer of the TOF and a corresponding signal 237 

is searched for. The channel with the center closest to the 23s 

track extrapolation point is selected if the distance is less 239 

than 10 cm. This rather weak criterion results in a high 240 

matching efficiency while keeping the fraction of wrongly 241 

associated tracks below 1% in the low-density environment 242 

presented by pp collisions. 243 

The dE/dx measurements in the ITS are used to iden- 244 

tify hadrons in two independent analyses, based on dif- 245 

ferent tracking algorithms. One analysis uses the ITS- 246 

TPC combined tracking, while the other is based on ITS 247 

stand-alone tracks. The combined ITS-TPC tracking re- 24s 

suit serves as a cross-check of both the ITS stand-alone 249 

and the TPC results in the overlap region. The ITS stand- 250 

alone analysis extends the acceptance to lower p t than the 251 
TPC or ITS-TPC analyses. 

The combined ITS-TPC analysis uses the same track 253 

selection criteria as the TPC only analysis, with the ad- 254 

ditional requirement of at least four clusters in the ITS, 255 

out of which at least one must be in the SPD and at least 256 

three in SSD+SDD. This further reduces the contamina- 257 

tion of secondaries and provides high resolution on track 25s 

impact parameter and optimal resolution on the dE/dx. 259 

The ITS stand-alone tracking uses a similar selection, with 260 

a different \ 2 selection and a different DCA selection. In 201 

the current tracking algorithm, ITS clusters are assigned 262 

a larger position error to account for residual misalign- 263 

ment of the detector. As a result, the x 2 values are not 264 

properly normalized, but the selection was adjusted to be 265 

equivalent to the TPC \ 2 selection by inspecting the dis- 266 

tributions. The DCA cut in the ITS analysis uses the same 26? 

Pt-dependent parametrization as for TPC tracks, but with 268 

different parameters to account for the different resolution. 259 
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270 2.4 Monte-Carlo Calculations 



271 The efficiency and other correction factors including ac- 

272 ceptance (jointly called efficiency in the following discus- 

273 sion) used in this paper are calculated from a Montc- 

274 Carlo simulation, based on over two million events pro- 

275 duced with the PYTHIA 6.4 event generator [TS] (tune 

276 D6T [16]), propagated through the detector with the GE- 

277 ANT3 [T7] transport code. Dead and noisy channels as 

278 well as beam position and spread have been taken into 

279 account. A simulation based on the PHOJET event gen- 

280 erator [TS] is also used as a cross check. 

281 GEANT3 is known to reproduce the absorption cross 

282 sections of hadrons incorrectly. The transport code FLU- 

283 KA contains a more accurate description of these cross 

284 sections [TM2"T] . and a dedicated simulation is used to 

285 calculate a correction to the GEANT3 efficiency calcula- 

286 tion [22) . This is relevant mainly for antiprotons at low 

287 p t , where the correction is on the order of 10%. For other 

288 particles and at higher p t , the difference between GEANT 

289 and FLUKA calculations is negligible. 



290 2.5 Particle Identification 



2.5.1 Particle identification in the ITS 322 

In both the ITS stand-alone and in the ITS-TPC analy- 323 

ses, the dE/dx measurement from the SDD and the SSD 324 

is used to identify particles. The stand-alone tracking re- 325 

suit extends the momentum range to lower p t than can be 326 

measured in the TPC, while the combined tracking pro- 327 

vides a better momentum resolution. 32s 

The energy loss measurement in each layer of the ITS 329 

is corrected for the track length in the sensitive volume 330 

using tracking information. In the case of SDD clusters, a 331 

linear correction for the dependence of the reconstructed 332 

raw charge as a function of drift time due to the com- 333 

bined effect of charge diffusion and zero suppression is 334 

also applied [5]. For each track, dE/dx is calculated using 335 

a truncated mean: the average of the lowest two points 336 

in case four points are measured, or a weighted sum of 337 

the lowest (weight 1) and the second lowest point (weight 33s 

1/2), in case only three points are measured. 339 




momentum [GeV/c] 



291 The dE/dx and TOF signals are used for particle iden- 

292 tification as a function of the momentum p, whereas the 

293 final spectra are given as a function of the transverse mo 

294 mentum p t . 

295 In the case of the TPC and ITS analyses, particles were 

296 identified via the specific energy loss dE/dx. Unique iden- 

297 tification on a track-by-track basis is possible in regions 

298 of momentum where the bands are clearly separated from 

299 each other. In overlapping areas, particle identification is 

300 still possible on a statistical basis using fits to the energy 

301 loss distribution in each p t -bin. The fits are performed on 

302 the distribution of the difference between the measured 

303 and the expected energy deposition for tracks within the 

304 selected rapidity range \y\ < 0.5. This compensates for 

305 the very steep slope of the Bethe-Bloch in the 1//3 2 region 

306 which would make the d-E/dx-distribution in a simple p t 

307 or p-slice non-Gaussian. The calculated expected energy 

308 loss depends on the measured track momentum p and the 

309 assumed mass for the particle. The procedure is therefore 

310 repeated three times for the entire set of tracks, assuming 

311 the pion, kaon, and proton mass. 

312 In the TPC analysis, the difference 

313 [dE/dx] mcas ~[dE/dx(pid,ptot)]caic is used. For the ITS 

314 the difference of the logarithm of the measured and calcu- 

315 lated energy deposit ln[cLE /dx mcas ]-ln[di? / dx{pid,p to t) ca \c] 

316 is taken to suppress the non-gaussian tails originating from 

317 the smaller number of dE/dx measurements. 

318 In the case of the TOF, the identification is based on 

319 the time-of-flight information. The procedure for the ex- 

320 traction of the raw yields differs slightly from the one used 

321 for TPC and ITS, and is described in Section [2.5. 31 



Fig. 1. (Color online) Specific energy loss dE/dx vs. momen- 
tum for tracks measured with the ITS. The solid lines are a 
parametrization (from [23]) of the detector response based on 
the Bethe-Bloch formula. 

Figure [T] shows the truncated mean dE/dx for the sam- 340 

pie of ITS stand-alone tracks along with the PHOBOS 341 

parametrization of the most probable value |23] . 342 

For the ITS stand-alone track sample, the histograms 343 

are fitted with three Gaussians and the integral of the 344 

Gaussian centered at zero is used as the raw yield of the 345 

corresponding hadron species. In a first step, the peak 346 

widths (j of the peaks are extracted as a function of pt 347 

for pions and protons in the region where their dE/dx 34s 

distributions do not overlap with the kaon (and electron) 349 

distribution. For kaons, the same procedure is used at low 350 

p t , where they are well separated. The p t -dependence of 351 

the peak width is then extrapolated to higher p t with the 352 

same functional form used to describe the pions and pro- 353 

tons. The resulting parametrizations of the pt dependence 354 

of a are used to constrain the fits of the In [dE/dx] distri- 355 

butions to extract the raw yields. 356 

For the ITS-TPC combined track sample, a non-Gau- 357 

ssian tail is visible. This tail is a remnant of the tail of the 35s 
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Fig. 2. (Color online) Distribution of ln[di5/da;] moas — ln[d_B/dx(K)] ca i c measured with the ITS in the two p t -ranges, 300-350 
MeV/c (upper panels) and 400-450 MeV/c (lower panels), using the kaon mass hypothesis. The left panels show the result for 
ITS-TPC combined tracks, while the right panels show the ITS stand-alone result. The lines indicate fits as described in the 
text. 



359 Landau distribution for energy loss. It was verified using 

360 simulations that the shape and size of the tail are com- 

361 patible with the expectations for a truncated mean using 

362 two out of four samples. The tail is not as pronounced 

363 for the ITS stand-alone track sample, due to the limited 

364 momentum resolution. The distribution is fitted with a 

365 combination of a Gaussian and an exponential function 

366 for the main peak and another exponential function to 

367 describe the tail of a background peak. This functional 

368 form provides an accurate description of the peak shape 

369 in the detector simulation, as well as the measured shape. 

370 Examples of dE/dx distributions are shown in Fig. [5] 

371 for negative tracks using the kaon mass hypothesis in two 

372 different pt intervals for both ITS stand-alone tracks (right 

373 panels) and ITS-TPC combined tracks (left panels). 

374 Efficiency correction The raw hadron yields extracted 

375 from the fits to the dE/dx distributions are corrected 

376 for the reconstruction efficiency determined from Monte- 

377 Carlo simulations, applying the same analysis criteria to 

378 the simulated events as to the data. Secondary particles 

379 from interactions in the detector material and strange par- 

380 tide decays have been subtracted from the yield of both 

381 simulated and real data. The fraction of secondaries af- 

382 ter applying the track impact-parameter cut depends on 

383 the hadron species and amounts to 1-3% for pions and 

384 5-10% for protons depending on p t . The secondary-to- 

385 primary ratio has been estimated by fitting the measured 

386 track impact-parameter distributions with three compo- 



nents, prompt particles, secondaries from strange particle 387 

decays and secondaries produced in the detector material 388 

for each hadron species. Alternatively, the contamination 389 

from secondaries have been determined using Monte-Carlo 390 

samples, after rescaling the A yield to the measured val- 391 

ues |24| . The difference between these two procedures is 392 

about 3% for protons and is negligible for other particles. 393 

Figure [3] shows the total reconstruction efficiency for 394 

primary tracks in the ITS stand-alone, including the ef- 395 

fects of detector and tracking efficiency, the track selection 396 

cuts and residual contamination in the fitting procedure, 397 

as determined from the Monte-Carlo simulation. This ef- 39s 

ficiency is used to correct the measured raw yields after 399 

subtraction of the contributions from secondary hadrons. 400 

The measured spectra are corrected for the efficiency of 401 

the primary vertex reconstruction with the SPD using 402 

the ratio between generated primary spectra in simulated 403 

events with a reconstructed vertex and events passing the 404 

trigger conditions. 405 

Systematic errors are summarized in Table [1] The sys- 406 

tematic uncertainty from secondary contamination has been 407 

estimated by repeating the full analysis chain with differ- 4os 

ent cuts on the track impact parameter and by comparing 409 

the two alternative estimates outlined above. The effect of 410 

the uncertainty in the material budget has been estimated 411 

by modifying the material budget in the Monte-Carlo sim- 412 

ulations by ±7%, which is the present uncertainty of the 413 

ITS material budget. The systematic contribution from m 

the fitting procedure to the ln[di?/d.x] mcas — ln[di?/da;(i)] ca i c 415 

distributions has been estimated by varying the fit condi- 416 
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Fig. 3. (Color online) Efficiency for pions, kaons and protons 
for the ITS stand-alone analysis as obtained from Monte-Carlo 
simulations. 



Table 1. Summary of systematic errors in the efficiency cor- 
rection of the ITS analysis. 



systematic errors 


7T± 


K± 


p and p 


secondary contamination 
from material 


negl. 


negl. 


negl. 


secondary contamination 
from weak decay 


< 1% 


negl. 


3% 


material budget 
highest p t bin 
lowest p t bin 


< 1% 
5% 


< 1% 
2% 


1% 
3% 


ITS efficiency 
all p t bins 
lowest p t bin 


2% 
12% 


2% 
13% 


2% 
11% 


ln(dS/da;) distr. 
fitting procedure 


1% 


5% 


3.5% 



tion and by comparing to an independent analysis using 
a track-by-track identification approach based on the dis- 
tance between the measured and expected dE/dx values 
normalized to its resolution. The residual imperfections in 
the description of the ITS detector modules and dead ar- 
eas in the simulation introduce another uncertainty in the 
ITS tracking efficiency. This is estimated by varying the 
cuts on the number of clusters and on the track \ 2 both 
in data and in Monte-Carlo simulations. 

In the lowest p t -bins, a larger systematic error has been 
assigned to account for the steep slope of the tracking effi- 
ciency as a function of the particle transverse momentum 
(see Fig. El). 



2.5.2 Particle identification in the TPC 

Particle identification is based on the specific energy de- 
posit of each particle in the drift gas of the TPC, shown in 
Fig. 2] as a function of momentum separately for positive 
and negative charges. The solid curves show the calibra- 
tion curves obtained by fitting the ALEPH parametriza- 
tion of the Bethe-Bloch curve [25 to the data points in 
regions of clear separation. 

The calibration parameters have mostly been deter- 
mined and tested via the analysis of cosmic rays. The 
pad-gain factors have been measured using the decay of 
radioactive fgKr gas released into the TPC volume (for a 
detailed description see [7]). 
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Fig. 4. (Color online) Specific energy loss dE/dx vs. momen- 
tum for tracks measured with the ALICE TPC. The solid lines 
are a parametrization of the Bethe-Bloch curve |25| . 



As in the case of the ITS, a truncated-mean procedure 
is used to determine dE/dx (60% of the points are kept). 
This reduces the Landau tail of the dE/dx distribution to 
the extent that it is very close to a Gaussian distribution. 

Examples of the dE/dx distribution in some pt bins 
are shown in Fig. [5] The peak centered at zero is from 
kaons and the other peaks are from other particle species. 
As the background in all momentum bins is negligible, the 
integrals of the Gaussian give the raw yields. 



Efficiency correction The raw hadron spectra are cor- 
rected for the reconstruction efficiency, shown in Fig. [6j 
determined by doing the same analysis on Monte-Carlo 
events. The efficiency is calculated by comparing the num- 
ber of reconstructed particles to the number of charged 
primary particles from PYTHIA in the chosen rapidity 
range. For transverse momenta above 800 MeV/c the ef- 
ficiency saturates at roughly 80%. For kaons, the decay 
reduces the efficiency by about 30% at 250 MeV/c and 
12% at 1.5 GeV/c. The range with a reconstruction effi- 
ciency lower than 60% (for pions and protons) is omitted 
for the analysis corresponding to a low-pt cut-off of 200 
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Fig. 6. (Color online) Efficiency of charged pions, kaons, and 
protons for the spectra extracted with the TPC. 
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Fig. 5. (Color online) Distribution of 

( [dE/ dx] mC as - [dE/ dx (kaon)] ca ic)/ [dS/da;(kaon)] ca i c mea- 
sured with the TPC for several p t -bins showing the separation 
power. The solid lines are Gaussian fits to the distributions. 



464 MeV/c for pions, 250 MeV/c for kaons, and 400 MeV/c 

465 for protons. 

466 Protons are corrected for the contamination of sec- 

467 ondaries from material and of feed down from weak de- 

468 cays. The feed down was determined by two independent 

469 methods. Firstly, the contamination obtained from Montc- 

470 Carlo simulation was scaled such that it corresponds to 

471 the measured yield of As in the data [H]. Secondly, the 

472 shape of the impact parameter distribution was compared 

473 to the Monte-Carlo simulation. Weak decays produce a 

474 non-Gaussian tail in the distribution of primary particles 

475 whereas secondaries from material generate a flat back- 

476 ground .22 . The remaining difference between the meth- 

477 ods is included in the systematic error. The correction for 



weak decays amounts to up to 14% and the correction for 478 

secondaries from material up to 4% for protons with 400 479 

MeV/c < pt < 600 MeV/c. For other particle species and 4so 

other transverse momenta the contamination is negligible. 48i 

The systematic errors in the track reconstruction and 482 

in the removal of secondary particles have been estimated 483 

by varying the number of standard deviations in the dis- 484 

tance-to- vertex cut, using a fixed cut of 3 cm instead of 485 

the variable one, and varying the SPD-TPC matching cut. 486 

Their impact on the corrected spectra is less than 5%. The 487 

influence of the uncertainty in the material budget has 488 

been examined by varying it by 7%. This resulted in the 489 

systematic errors given in Table [2] The uncertainty due 490 

to a possible deviation from a Gaussian shape has been 491 

established by comparing the multi-Gauss fit with a 3-cr 492 

band in well separated regions. The precision of the kink 493 

rejection is estimated to be within 3%. 494 

The correction for the event selection bias has been 495 

tested with two event generators, PYTHIA [TSIIH] and 496 

PHOJET [TH] and the corresponding uncertainty is less 497 

than 1%. 498 



2.5.3 Particle identification with the TOF 499 

Particles reaching the TOF system are identified by mea- 500 

suring their momentum and velocity simultaneously. 501 

The velocity j3 = L/ixoF is obtained from the mea- 502 

sured time of flight tjtof an d t ne reconstructed flight path 503 

L along the track trajectory between the point of closest 504 

approach to the event vertex and the TOF sensitive sur- 505 

face. The measured velocities are shown as a function of 506 
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Table 2. Summary of systematic errors in the efficiency cor- 
rection in the TPC analysis. 



systematic errors 




K± 


p and p 


secondary contamination 
from material 


negl. 


negl. 


< z/o 


secondary contamination 
from weak decay 


< 4% 




< 10% 


energy loss and 
absorption in material 


< 1% 


< 1% 


< 2% 


kink rejection 


negl. 


< 3% 




non-Gaussianity of 
dE/dx signal 


negl. 


negl. 


negl. 


matching to ITS 




< 3% 
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Fig. 7. (Color online) j3 of tracks of particles measured by 
TOF vs. their momentum. 



507 the momentum p at the vertex in Fig. [7] The bands corre- 

508 sponding to charged pions, kaons and protons are clearly 

509 visible. The width of the bands reflects the overall time- 

510 of-fiight resolution of about 180 ps, which depends on the 

511 TOF timing signal resolution, the accuracy of the recon- 

512 structed flight path and the uncertainty of the event start 

513 time, tff. This last contribution is related to the uncer- 

514 tainty in establishing the absolute time of the collision. 

515 In the present sample this fluctuated with respect to the 

516 nominal time signal from the LHC with a a of about 140 

517 ps due to the finite size of the bunches. 

sis To improve the overall time-of-flight resolution, the 

519 TOF information itself is used to determine ig" in events 

520 having at least three tracks with an associated TOF sig- 

521 nal. This is done with a combinatorial algorithm which 

522 compares the TOF times with the calculated times of the 

523 tracks for each event for different mass hypotheses. Us- 

524 ing this procedure, the start-time has been improved for 

525 44% of the tracks having an associated TOF signal and 

526 is rather independent on the momentum of the tracks. In 

527 this way the precision on the event start-time is about 85 

528 ps on average. 



Finally, tracks whose particle identity as determined 529 

from the TOF information is not compatible with the one 530 

inferred from the dE/dx signal in the TPC within five a 531 

have been removed. This TOF-TPC compatibility crite- 532 

rion rejects about 0.6% of the tracks and further reduces 533 

the small contamination coming from tracks incorrectly 534 

associated with a TOF signal. 535 
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Fig. 8. (Color online) Distribution of the time difference be- 
tween the measured TOF signal and the average of the cal- 
culated times for pions and kaons for several p t -bins for posi- 
tively charged particles. The fits are performed using Gaussian 
shapes. 
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For each particle species i, the expected time of flight 
i* alc is calculated by summing up the time-of-flight incre- 
ments Atk — Alk-\/pi + Tn3/pk at each tracking step, with 
Pk being the local value of the track momentum, the 
mass of the particle, and Alk the track-length increment 
along its trajectory. The yields of 7r, K and p are obtained 
from the simultaneous fit of the distribution of the time 
difference S between measured £tof and the average be- 
tween the calculated time for pions and kaons 



S = (t W + t K )calc/2 — ^TOF- 



(1) 



The symmetric treatment of kaons and pions in the defi- 
nition of 5 ensures that the kaon and pion peak are both 
Gaussian. Extracting the yield for different species in a si- 
multaneous fit guarantees that the resulting number of pi- 
ons, kaons and protons matches the total number of tracks 
in the given momentum bin. 

The distribution of the variable S is shown in Fig. [5] for 
three different transverse momentum bins for positive par- 
ticles. The curves show the results of the three- Gaussian 
fit used to extract the raw yields. The integral of the fit 
result has been constrained to the number of entries in the 
distribution, and the means and the widths are allowed to 
vary within 5% and 10%, respectively, of their nominal 
values. The only free parameters in the fit are therefore 
the relative normalizations between the Gaussians. 

The raw yields are extracted in different p t -bins using 
a rapidity selection |y p | < 0.5, where y p is the rapidity 
calculated with the proton mass. For pions and kaons, 
this condition results in a larger y-acceptance and in both 
cases, the fraction outside of \y\ < 0.5 has been subtracted 
in each p t -bin taking into account the y-distribution of the 
yields within the pions and kaons peaks. 



Efficiency correction Since the track selection used in the 
TOF analysis is the same as the one described in the 
TPC analysis fsubsection 12.5.2]) . the same tracking and 
feed-down corrections are applied. In the case of the TOF 
analysis, an additional correction is needed in order to take 
into account the fraction of the particles reconstructed by 
the TPC with an associated signal in TOF. This matching 
efficiency includes all sources of track losses in the propa- 
gation from the TPC to the TOF (geometry, decays and 
interactions with the material) and its matching with a 
TOF signal (the TOF intrinsic detector efficiency, the ef- 
fect of dead channels and the efficiency of the track-TOF 
signal matching procedure). The TOF matching efficiency 
has been derived from Monte-Carlo events as the fraction 
of TPC reconstructed tracks having an associated TOF 
signal and is shown in Fig. [5] for each hadron species. The 
main factors limiting the TOF matching efficiency are the 
loss due to geometrical acceptance (w 15%), the number 
of dead or noisy channels 10%) and the absorption of 
particles in the material of the transition radiation detec- 
tor (« 8%). 

The TOF matching efficiency has been tested with 
data, using dE/dx in the TPC to identify the particles. 
Good agreement between the efficiencies obtained from 
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Fig. 9. (Color online) The TOF matching efficiency is shown 
for the three particles, separately, for (top) positive and (bot- 
tom) negative particles. 

Table 3. Summary of systematic errors in the TOF analysis. 



systematic errors 




K± 


p and p 


TOF 


< 3% 


< 6% 


< 4% 


matching 






(pt> 1 GeV/c) 


efficiency 






< 7.5% 








(pt= 0.7 GeV/c) 


PID procedure 


< 2% 


< 7% 


< 3% 



the data and from Monte-Carlo simulations is observed in 
case of pions and kaons, with deviations at the level of, 
at most, 3% and 6% respectively, over the full transverse- 
momentum range. The observed differences are assigned 
as systematic errors, see Table |31 In the case of protons 
and antiprotons, larger differences are observed at pt be- 
low 0.7 GeV/c, where the efficiency varies very rapidly 
with momentum. This region is therefore not considered 
in the final results (see Table [3]) ■ 

Other sources of systematic errors related to the TOF 
PID procedure have been estimated from Monte-Carlo 
simulations and cross-checked with data. They include the 
effect of the residual contribution from tracks wrongly as- 
sociated with TOF signals, and the quality and stability 
of the fit procedure used for extracting the yields. Tabled 
summarizes the maximal value of the systematic errors ob- 
served over the full transverse momentum range relevant 
in the analysis, for each of the sources mentioned above. 



2.6 Kaon Identification using their decay within the 
TPC 

In this section, the determination of the yields of charged 
kaons identified by their weak decay (kink topology) inside 
the TPC detector is described. These tracks are rejected 
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in the previously described TPC analysis. This procedure 
allows an extension of the study of kaons to intermedi- 
ate momenta, on a track-by-track level, although in this 
analysis the pt reach is limited by statistics. 

The kinematics of the kink topology, measured as a 
secondary vertex with one mother and one daughter track 
of the same charge, allows the separation of kaon decays 
from the main source of background kinks coming from 
charged pion decays. The decay channels with the highest 
branching ratio (B.R.) for kaons are the two-body decays 



(1) K± 

(2) K± 



, (B.R. 63.55%) 
, (B.R. 20.66%) 



Three-body decays with one charged daughter track 
(B.R. 9.87%) as well as three-body decays into three char- 
ged pions (B.R. 5.6%) are also detected. 

The algorithm for reconstructing kinks as secondary 
vertices is applied inside a fiducial volume of the TPC with 
radius 120 cm < R < 210 cm in order to have a minimum 
number of clusters for reconstructing both the mother and 
daughter tracks. Inside this volume a sufficient number of 
kinks can be found since the cr of kaon and pion decays 
are 3.7 m and 7.8 m, respectively. The mother track of 
the kink has been selected with similar criteria to those 
of the TPC tracks used for the dE/dx analysis, except 
that the minimum required number of clusters per track 
is 30, because the kink mother track does not traverse the 
entire TPC. The relation between the number of clusters 
per mother track and the radius R of the kink is used as 
a quality check of the kink reconstruction procedure. 

The identification of kaons from kink topology and its 
separation from pion decay is based on the decay kine- 
matics. The transverse momentum of the daughter with 
respect to the mother's direction, q t , has an upper limit 
of 236 MeV/c for kaons and 30 MeV/c for pions for the 
two-body decay to fi + z^. The corresponding upper limit 
for the two-body decay (2) K -> tt + tt° is 205 MeV/c. 
All three limits can be seen as peaks in Fig. [TU] (a), which 
shows the qt distribution of all measured kinks inside the 
selected volume and rapidity range \y\ < 0.7. Selecting 
kinks with q t > 40 MeV/c removes the majority of ir- 
decays as shown by the dashed (before) and solid (after) 
histograms. 

The invariant mass for the decay into /i^ + is cal- 
culated from the measured difference between the mother 
and daughter momentum, their decay angle, assuming zero 
mass for the neutrino. Figure [10] (b) shows the invariant 
mass for the full sample of kinks (dashed line) and for 
the sample after applying the preceding cuts (full line). 
The masses of pions and kaons are reconstructed at their 
nominal values. The third peak at 0.43 GeV/c originates 
from the K — > ir + ir° decay for which the invariant mass 
is calculated with wrong mass assumptions for the daugh- 
ter tracks. The broad structure originates from three-body 
decays of kaons. 

At this stage, we have a rather clean sample of kaons 
as demonstrated in Fig. [10] (c) showing the dE/dx vs. the 
mother momentum. Most of the tracks are within a 3.5er 
band with respect to the corresponding Bethe-Bloch curve 
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Fig. 10. (Color online) (a) q t distribution of the daugh- 
ter tracks with respect to mother momentum for all recon- 
structed kinks inside the analyzed sample. The dashed(solid) 
histograms show the distribution before (after) applying the 
qt > 40 MeV/c cut. (b) Invariant mass of the two-body decays 
K ± /-k ± — > /i* + U/j, for candidate kaon kinks. Solid curve: after 
applying q t >40 MeV/c; dashed curve: without this selection 
(hence also showing the pion decays), (c) dE/dx of kinks as a 
function of the mother momentum, after applying the full list 
of selection criteria for their identification. 



of kaons. The few tracks outside these limits are at mo- 
menta below 600 MeV/c (less than 5%) and they have 
been removed in the last analysis step. 



Efficiency and acceptance The total correction factor in- 
cludes both the acceptance of kinks and their efficiency 
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Fig. 11. (Color online) Upper panel: The acceptance of kaons 
decaying in the fiducial volume of the TPC as a function of the 
kaon p t for K + (full-triangles) and K~ (open-squares). Lower 
panel: The efficiency of reconstructed kaons from kinks as a 
function of the p t (mother), separately for K + (full-triangles) 
and K~ (open-squares). The contamination from wrongly as- 
sociated kinks is also plotted for both charges (lower set of 
points). 



677 (reconstruction and identification). The study has been 

678 performed for the rapidity interval \y\ < 0.7, larger than 

679 the corresponding rapidity interval for the other studies 

680 in order to reduce the statistical errors. 

681 The acceptance is defined as the ratio of weak decays 

682 (two- and three-body decays) whose daughters are inside 

683 the fiducial volume of the TPC to all kaons inside the same 

684 rapidity window (Fig. [TTJ upper part). It essentially re- 
ess fleets the decay probability. However, the acceptance is not 

686 the same in the low-momentum region for both charges of 

687 kaons, since the interaction cross section of the negative 
ess kaons with the ITS material is higher than that of the pos- 

689 itive kaons. As a result, the acceptance of positive kaons 

690 is larger at low momenta. 

691 The efficiency is the ratio of reconstructed and identi- 

692 fied kaons divided by the number of kaon decays within 

693 the acceptance as shown in Fig. 1111 (lower part), as a func- 

694 tion of the kaon p t . It reaches about 60% at 0.7 GeV/c 

695 and decreases gradually at higher transverse momenta, as 

696 the angle between mother and daughter tracks becomes 

697 smaller. The decay angle of kaon kinks allows their iden- 



tification up to high momenta, e.g. at p t of 5 GeV/c the 
values are between 2° and 15°. 

The contamination due to random associations of pri- 
mary and secondary charged tracks has been established 
using Monte-Carlo simulations and it is systematically 
smaller than 5% in the studied pt-range as also shown 
in Fig. 1111 Hadronic interactions are the main source of 
these fake kinks (65%). 

The systematic error due to the uncertainty in the ma- 
terial budget is about 1% as for the TPC analysis. The 
quality cuts remove about 8% of all real kaon kinks, which 
leads to a systematic error of less than 1%. The main un- 
certainty originates from the efficiency of the kink finding 
algorithm which has an uncertainty of 5%. 



3 Results 

Figure [T2"l shows a comparison between the results from the 
different analyses. The spectra are normalized to inelastic 
collisions, as explained in Sec. 12.21 The kaon spectra ob- 
tained with various techniques, including spectra |24j . 
are compared in Fig. 1131 The very good agreement demon- 
strates that all the relevant efficiencies are well reproduced 
by the detector simulation. 

The spectra from ITS stand-alone, TPC and TOF are 
combined in order to cover the full momentum range. The 
analyses from the different detectors use a slightly differ- 
ent sample of tracks and have largely independent sys- 
tematics (mainly coming from the PID method and the 
contamination from secondaries). The spectra have been 
averaged, using the systematic errors as weights. From this 
weighted average, the combined, pt-dependent, systematic 
error is derived. The combined spectra have an additional 
overall normalization error, coming primarily from the un- 
certainty on the material budget (3%, Sec. I2.5|) and from 
the normalization procedure (2%, Sec. 12. 2p . 

The combined spectra shown in Fig. [T?] are fitted with 
the Levy (or Tsallis) function (see e.g. [26ll27| ) 



d 2 7V 
dp t dy 



dN 



Ptx 



(n- l)(n-2) 



dy nC(nC + mo(n — 2)) 



m t - rap 
nC 



(2) 

with the fit parameters C, n and the yield dN /dy. This 
function gives a good description of the spectra and has 
been used to extract the total yields and the (pt), summa- 
rized in Table |3J The x 2 /degree-of- freedom is calculated 
using the total error. Due to residual correlations in the 
point-by-point systematic error, the values are less than 1. 
Also listed are the lowest measured pt-bin and the fraction 
of the yield contained in the extrapolation of the spectra to 
zero momentum. The extrapolation to infinite momentum 
gives a negligible contribution. The systematic errors take 
into account the contributions from the individual detec- 
tors, propagated to the combined spectra, the overall nor- 
malization error and the uncertainty in the extrapolation. 
The latter is evaluated using different fit functions (mod- 
ified Hagedorn [35] and the UA1 parametrization or 
using a Monte-Carlo generator, matched to the data for 
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Fig. 12. (Color online) Transverse momentum spectra 
d 2 N / (dptdy) for \y\ < 0.5 of positive (upper part) and negative 
(lower part) hadrons from the various analyses. Only system- 
atic errors are plotted. 



Pt < 1 GeV/c (PYTHIA [IS], with tunes D6T [IB], CSC 
and PerugiaO [3D], or PHOJET [IB]). While none of these 
alternative extrapolations provides a description as good 
as the one from the Levy fit, we estimate from this pro- 
cedure an uncertainty of about 25% of the extrapolated 
part of the yield. 

The ratios of 7r + /vT~and K + /K~ as a function of pt are 
close to unity within the errors, allowing the combination 
of both spectra in the Levy fits. The p/p ratio as a function 
of pt has been studied with high precision in our previous 
publication (22| . It is pt- independent with a mean value of 
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Fig. 13. (Color online) Comparison of charged kaon spectra, 
obtained from the combined ITS stand-alone, TPC, TOF anal- 
ysis, from the kink topology and K° spectra from Ref. [24] 
. Only statistical errors are shown. 



0.957±0.006(stat)±0.014(syst). Also here we used the sum to 

of both charges. Table [5] summarizes the fit parameters 762 

along with the yields and mean p t . The errors have been 763 

determined as for the individual fits. 764 

Our values on yield and (pt) given in Table [D and [SJ 765 

agree well with the results from pp collisions at the same 766 

y/s [3 1J . Figure [TS1 compares the (pt) with measurements 767 

in pp collisions at y/s = 200 GeV [32][33] and in pp re- 768 

actions at y/s = 900 GeV [31]. The mean p t rises very ?69 

little with increasing ^/s despite the fact that the spectral 770 

shape clearly shows an increasing contribution from hard 771 

processes. It was already observed at RHIC that the in- 772 

crease in mean p t at y/s= 200 GeV compared to studies at 773 

y/s= 25 GeV is small. The values obtained in pp collisions 774 

are lower than those for central Au+Au reactions at \fs= 775 
200 GeV [3J. 

The spectra presented in this paper are normalized 777 

to inelastic events. In a similar study by the STAR Col- 77s 

laboration the yields have been normalized to NSD colli- 779 

sions [32 . In order to compare these two results, the yields 7so 

in Table[?]have been scaled to NSD events, multiplying by 7m 

1.185 (see Section 2.2). The yields of pions increase from 782 

Vs= 200 GeV to 900 GeV by 23%, while K+ rises by 45% 7 83 

and K- by 48%. 7 84 

Figure [TB] shows the K/7T ratio as a function of y/s both 785 

in pp (full symbols, [32 ( 34, 35 ) and in pp (open symbols, 786 

36 38.) collisions. For most energies, (K + +K _ ) /(7r + +7r~) 78? 

is plotted, but for some cases only neutral mesons were 788 

measured and K°/7r° is used instead. The pt-integrated im 

(K + +K _ )/(7r + +?T _ ) ratio shows a slight increase from 790 

Vs= 200 GeV (K/tt = 0.103 ± 0.008) to Vs= 900 GeV 791 

(K/tt^O.123 ± 0.004 ± 0.010) 32j, yet consistent within 792 

the error bars. The results at 7 TeV will show whether 793 
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Table 4. Integrated yield AN /Ay (\y\ < 0.5) with statistical and systematic errors, and (pt), as obtained from the fit with the 
Levy function together with the lowest p t experimentally accessible, the fraction of extrapolated yield and the x 2 / n <Af of the fit 
(see text). The systematic error of AN /Ay and of the (pt) includes the contributions from the systematic errors of the individual 
detectors, from the choice of the functional form for extrapolation and from the absolute normalization. 



Particle 


AN /Ay 


( Pt ) (GeV/c) 


Lowest p t (GeV/c) 


Extrapolation 


X 2 /ndf 


TT + 


1.493 ±0.004 ±0.074 


0.404 ± 0.001± 0.02 


0.10 


10% 


14.23/30 


TT~ 


1.485 ±0.004 ±0.074 


0.404 ± 0.001± 0.02 


0.10 


10% 


12.46/30 


K+ 


0.183 ± 0.004 ±0.015 


0.658 ± 0.006± 0.05 


0.20 


13% 


12.71/24 


K~ 


0.182 ± 0.004 ±0.015 


0.642 ± 0.006± 0.05 


0.20 


13% 


6.23/24 


P 


0.083 ± 0.002 ± 0.006 


0.768 ± 0.008± 0.06 


0.35 


21% 


13.79/21 


P 


0.079 ±0.002 ±0.006 


0.760 ± 0.008± 0.06 


0.35 


21% 


13.46/21 



Table 5. Results of the Levy fits to combined positive and negative spectra. See text and the caption of Table [4] for details on 
the systematic errors. 

Particle d7V/dy C (GeV) n ( p t ) (GeV/c) X 2 /ndf 

7T+ + tt- 2.977 ±0.007 ±0.15 0.126 ± 0.0005 ± 0.001 7.82 ± 0.06 ± 0.1 0.404 ± 0.001 ± 0.02 19.69/30 
K+ ± K~ 0.366 ± 0.006 ± 0.03 0.160 ± 0.003 ± 0.005 6.08 ±0.2 ±0.4 0.651 ± 0.004 ± 0.05 8.46/24 
p±p 0.162 ±0.003 ± 0.012 0.184 ± 0.005 ± 0.007 7.5 ±0.7 ±0.9 0.764 ± 0.005 ± 0.07 15.70/21 



794 the K/-7T ratio keeps rising slowly as a function of ^fs or 

795 saturates. 

796 Protons and antiprotons in Table [4] have been cor- 

797 rected for feed down (mainly from A), while the results 

798 from the STAR Collaboration are not. The proton spec- 

799 tra measured by PHENIX, on the other hand, have a lower 
soo pt-cut of 0.6 GeV/c. This makes a direct comparison with 

801 RHIC data difficult. 

802 Figure [T7] shows a comparison of the measured pion, 

803 kaon and proton spectra with several tunes of the PYTHIA 

804 event generator [15] and with PHO JET [JJ] . The PYTHIA 
80s CSC 306 tune provides a very poor description of the 

806 particle spectra for all species. Similar deviations were 

807 already seen for the unidentified charged hadron spec- 

808 tra [it]. The other PYTHIA tunes, PerugiaO [30] and 

809 D6T 16J, and PHO JET give a reasonable description of 
sio the charged pion spectra, but show large deviations in the 
an kaon and proton spectra. The measured kaon p t -spectrum 

812 falls more slowly with increasing p t than the event genera- 

813 tors predict. A similar trend is seen for the proton spectra, 

814 except for PYTHIA tune D6T, which describes the proton 
sis spectra reasonably well. 

816 The upper panel of Figure [T51 shows the pt-dependence 

817 of the K/7T and also the measurements by the E735 [35] 
sis and STAR Collaborations 32 . It can be seen that the 
aw observed increase of K/7r with p t does not depend strongly 

820 on collision energy. 

821 A comparison with event generators shows that at pt> 

822 1.2 GeV/c, the measured K/7T ratio is larger than any of 

823 the model predictions. It is interesting to note that while 

824 the spectra in the CSC tune are much steeper than the 

825 other tunes, the pt-dependence of the K/tt ratio is very 

826 similar. In the models, the amount of strangeness produc- 

827 tion depends on the production ratios of gluons and the 

828 different quark flavours in the hard scattering and on the 

829 strangeness suppression in the string breaking. The latter 



could probably be tuned to better describe the data. A 830 

similar disagreement between measured strangeness pro- 83i 

duction and PYTHIA predictions was found at RHIC en- 832 

ergies gQ]. 833 

In the bottom panel of Figure 1181 the measured p/tt 

ratio is compared to results at \fs— 200 GeV from the 835 

PHENIX Collaboration [UJ . Both measurements are feed- 836 

down corrected. At low p t , there is no energy-dependence 837 

of the p/7r ratio visible, while at higher p t > 1 GeV/c, the 838 

p/ir ratio is larger at \/~s= 900 GeV than at \fs~ 200 GeV 839 

energy. 840 

Event generators seem to separate into two groups, 84i 

one with high P /tt ratio (PYTHIA CSC and D6T), which 84 2 

agree better with the data and one group with a lower 843 

p/ir ratio (PHO JET and PYTHIA PerugiaO), which are 844 

clearly below the measured values. These comparisons can 845 

be used for future tunes of baryon production in the event 846 

generators. en 



4 Summary 848 

We present the first analysis of transverse momentum spec- 849 

tra of identified hadrons, 7r + , 7t~, K + , K~, p, and p in pp sso 

collisions at ^/s = 900 GeV with the ALICE detector. The ssi 

identification has been performed using the dE/dx of the 852 

inner silicon tracker, the dE/dx in the gas of the TPC, 853 

the kink topology of the decaying kaons inside the TPC 854 

and the time-of-flight information from TOF. The combi- 855 

nation of these techniques allows us to cover a broad range sse 

of momentum. 857 

Agreement in the K/tt ratio is seen when comparing sss 

to pp collisions at the Tevatron and SppS. Comparing our 859 

results with a similar measurement from the STAR Col- seo 

laboration using pp collisions at y/s= 200 GeV the shape sei 

of the spectra shows an increase of the hard component, 862 
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Fig. 14. (Color online) Transverse momentum spectra of pos- 
itive (top) and negative (bottom) hadrons from pp collisions 
at yfs= 900 GeV. Grey bands: total p t -dependent error (sys- 
tematic plus statistical); normalization systematic error (3.1 
not shown. The curves represent fits using a Levy function. 



but we observe only a slight increase of the mean p t -values. 
Whether the fraction of strange to non-strange particles 
rises with increasing yfs remains open until data at 7 TeV 
become available. 
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Fig. 15. (Color online) Mean p t as a function of the mass of the 
emitted particle in pp collisions at 900 GeV (ALICE, red solid 
circles, statistical and systematic errors) compared to results 
at yfs = 200 GeV (star markers, average values of the results 
from the STAR and the PHENIX Collaborations [32][33]) and 
pp reactions at yfs = 900 GeV [31] (open squares). Some data 
points are displaced for clarity. 
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Fig. 16. (Color online) Ratios (K + +KT)/(7r + +7r-) and K° /vr 
as a function of yfs. Data (full symbols) are from pp collisions, 
(at yfs = 17.9 GeV by NA49 [SUES], at y/s = 200 GeV by 
STAR [32] and at y/s = 900 ALICE, present work) and (open 
symbols) from pp interaction (at y/s = 560 GeV by UA5 [37] 
and at the TEVATRON by E735 36, 3H]). 
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Fig. 17. (Color online) Comparison of measured pion, kaon 
and proton spectra at yfs = 900 GeV (both charges combined) 
with various tunes of event generators. Statistical errors only. 
See text for details. 
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Fig. 18. (Color online) Ratios of (K++ K")/(tt + + tt") (up- 
per panel) and (p + p) / (tt + + -k~) (lower panel) as a function 
of pt from pp collisions at ^J~s= 900 GeV (statistical errors 
only). Values from the E735 Collaboration [35] and the STAR 
Collaboration [35] (upper part) and from the PHENIX Collabo- 
ration [41] (lower part) also are given. The dashed and dotted 
curves refer to calculations using PYTHIA and PHOJET at 
^/s= 900 GeV. 



